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CN Abstract 

We study the production of two electro-weak vector bosons in association with a jet via gluon 
fusion. In particular we consider the production of 'fZg, ZZg and W + W~g at hadron colliders 
and compute their cross-sections. Such processes have already produced large number of events at 
D the Large Hadron Collider (LHC) by now. These processes can be background to the Higgs boson 

production and a number of beyond the standard model scenarios. Therefore it is important to know 
00 the values of their contribution. The calculation is based on conventional Feynman diagram approach. 

In particular we find that the process gg — >• ZZg can make a significant contribution to the process 
pp —t ZZ + jet. 

' 1 Introduction 

The search for new physics at the LHC is in progress and the collider is delivering data presently at 
8 TeV centre-of-mass energy. The discovery of a fundamental scalar particle (most probably a Higgs 
boson) of mass around 125 GeV has received a lot of world-wide attention [!]. We expect more good 
news from the experiments at the LHC before the collider goes for a two years long pause. So far, the 
standard model (SM) of particle physics seems to be in excellent agreement with the collected data (more 
Q\ than 10 fb -1 ). There has been searches for the hints of physics beyond the SM such as supersymmetry, 
large extra dimensions, etc. But, as of now, there is no clear evidence [2]. The process of identifying the 
discovered fundamental scalar particle as the Higgs boson is also continuing. With the lack of signals for 
^vq beyond the SM scenarios, there is a need to look for the SM processes that were not accessible earlier at 
the Tevatron. Most of such processes have several particles in the final state, and/or occur at the higher 
order. Such processes can also contribute to the background to the new physics signals. One such class 
of processes is multi vector boson production in association with one or more jets. 

At the LHC centre-of-mass energy, the collider has another useful feature. In the proton-proton 
collisions, the gluon luminosity can be quite significant. It can even dominate over the quark luminosity 
in certain kinematic domains. Therefore, at the LHC, loop mediated gluon fusion processes can be 
important. Di- vector boson production via gluon fusion have been studied by many authors [31 01 [5j 
El El El El dni HU Lia HH US H6]. We consider another class of processes gg -»■ VV'g, where V and V 
can be any allowed combination of electro- weak vector bosons. These processes can be a background to 
Higgs boson production as well as new physics scenarios such as technicolor. At the leading order, these 
processes receive contribution from quark loop diagrams. The prototype diagrams are displayed in Fig 
([!]) . The computation for the process gg — > 775 has already been performed |17| [TB] . Preliminary results 
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for gg —> jZg were presented in [19] . Recently, Melia et al. have presented calculation for gg — > W + W~g 
|20| . In this brief report we present our results for the production of "fZg, ZZg and W + W~g at hadron 
colliders. Detailed analysis of these processes will be reported in future publications. 




Figure 1: The prototype diagrams for the processes gg — > VV'g. The wavy lines represent the appropriate 
combination of the 7, W, or Z boson. The last two classes (c) and (d) are relevant to WW g production 
only. We do not consider diagrams involving Higgs boson for the ZZg and WW g cases. 



2 The Process: gW 

The process gg — > VV'g contribute to the cross-section of the process pp — > VV' + jet at NNLO in a s . 
As discussed, although they occur at the one-loop level, their cross-sections can be significant due to 
large gluon luminosity at the LHC. The processes jZg and ZZg receive contribution from two main 
classes of quark-loop diagrams - pentagon and box type, as shown in (a) and (b) of Fig Q. Box class 
of diagrams are due to the triple gluon vertices; they can be further divided into three sub-classes. This 
sub-classification has it's own physical importance. These are separately gauge invariant with respect 
to the electro-weak vector bosons. Other diagrams can be obtained by suitable permutation of external 
legs. For each quark flavor, there are 24 pentagon-type and 18 box-type diagrams. Due to Furry's 
theorem, only half of the 42 diagrams are independent. Also pentagon diagrams give both vector as 
well as axial-vector contributions while box diagrams give only vector contribution. We work with five 
massless quark flavors and expect decoupling of the top quark. 

In the case of WW g process, instead of single quark flavor, two quark flavors of a single generation 
contribute to the above discussed pentagon and box diagrams. To keep matter simple, for this process, 
we work with first two generations of massless quarks. It is expected that the contribution from the third 
generation would not be significant in low Pt region |16| . There are additional box and triangle class of 
diagrams due to j/Z W + W~ vertex for each quark flavor. These are shown in (c) and (d) of Fig. 0. Due 
to Furry's theorem, the triangle diagrams with ^W + W~ coupling don't contribute and only axial part 
of the triangle diagrams with ZW + W~ coupling contribute to the amplitude. Since the axial coupling of 
Z— boson to quark is proportional to Tg value, the axial contributions from additional triangle and box 
diagrams, when summed over massless quark generation, vanish. The vector contribution of additional 
box- type diagrams is separately gauge invariant. Because of its color structure, it interferes with the 
axial part of the pentagon amplitude. We have explicitly checked that its contribution towards total 
cross-section is very small; therefore we have dropped this contribution. Thus, effectively we are left 
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with ZZg like contributions for WW g. The Higgs boson interference effects for the cases of ZZg and 
WW g are ignored in the present calculation. Our one loop processes, being the leading order processes, 
are expected to be finite, i.e. free from UV and IR divergences. 

The amplitude for our process has the following general structure 



tabc jabc 

M abc (gg^VV'g) = i J —My{VV'g) + —M A {VV'g), 
M v {VV'g) = - e 2 gl C v (VV'g) (Vy - By) , 

M A (Wg) = - e 2 g 3 s C A (VV'g) (V A ) • (1) 

Mv,A are amplitudes for the vector and axial- vector parts of the full amplitude under consideration. 
Because of Bose symmetry M.y — > —My under exchange of any two external gluons while M A remains 
same. By and Vy a are box and pentagon contributions from a single flavor (single generation for WWg 
case) of quarks. The structure of the amplitude suggests that the vector and axial-vector contributions 
should be separately gauge invariant. Moreover due to the color structure when we square the amplitude, 
the interference between the axial and vector contribution vanishes, i.e. 

\M(gg VV'g)\ 2 = (g\Mv\ 2 + y \M A \ 2 ^j . (2) 

Therefore the cross-section of the process is an incoherent sum of the vector and axial- vector contri- 
butions. The couplings Cy tA for various cases are listed below. Contributions from all the relevant quark 
flavors (described above) are included appropriately. 
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C A (WWg) = *-(-).). (3) 
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3 Calculation and Numerical Checks 

For each class of diagrams, we write down the prototype amplitudes using the SM Feynman rules 
The amplitude of all other diagrams are generated by appropriately permuting the external momenta and 
polarizations in our code. The quark loop traces without 75 are calculated in n-dimensions while those 
with 75 are calculated in 4-dimensions using FORM [22]. We do not need any n-dimensional prescription 
for 75 as both pentagon and box diagrams are free of anomaly. The amplitude contains tensor loop 
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integrals. In the case of pentagon-type diagrams, the most complicated integral is rank-5 tensor integral 
(E^f" 75 ); while for the box-type dia grams, rank-4 tensor integral (D pvpu ) is the most complicated one 

„ & = f _d^k_ k»k»kPk°k s 

J (2vr)« A^iA^A^ ' 1 ' 

D „u P a = f d n k k»k"kPkT 

J (2tt)» N NiN 2 N 3 " 1 ' 

Here, N = kf — m q + ie and ki is the momentum of the i th internal line in the corresponding scalar 
integrals, n = (4 — 2e) and m q is the mass of the quark in the loop. Five point tensor and scalar 
integrals are written in terms of box tensor and scalar integrals using 4-dimensional Schouten Identity. 
For example, a five point scalar integral, in 4-dimensions, can be expressed in terms of five box scalar 
integrals [2"4] 

4 

£0(0,1,2,3,4) = £c i D< i) , (6) 

i=0 

where is the box scalar integral obtained after removing the i th propagator in Eq. Any 0(e) 
correction to the above relation in n = (4 — 2e)- dimensions does not contribute to our amplitude. 
This is because the contributions of the pentagon-type diagrams are UV finite and we have regulated 
mass singularities by giving small mass to the quarks. The box tensor integrals are reduced into the 
standard scalar integrals - Aq, Bq, Co and -Do using FORTRAN routines that follows from the reduction 
scheme developed by Oldenborgh and Vermaseren [23]. Thus we require box scalar integrals with two 
massive external legs, at the most. The scalar integrals with massive internal lines are calculated using 
OneLoop library [25J. Because of very large and complicated expression of the amplitude, we calculate 
the amplitude numerically before squaring it. This requires numerical evaluation of the polarization 
vectors of gauge bosons. We choose real basis, instead of helicity basis for the polarization vectors to 
calculate the amplitude. This is to reduce the size of compiled program and the time taken in running 
the code. We have used RAMBO to generate three particle phase-space for our processes [26]. 

The processes gg — > VV'g are leading order one- loop processes. The full amplitude in (1) should 
be both UV as well as IR finite. However, individual diagrams may be UV and/or IR divergent. IR 
divergence is relevant to only light quark cases. All these singularities are encoded in various scalar 
integrals. To make UV and IR finiteness check on our amplitude we have derived all the required scalar 
integrals (up to box scalar integrals with two massive external legs) analytically following t'Hooft and 
Veltman |27j . We regulate the UV divergence of the scalar integrals using dimensional regularization and 
infrared singularities by using small quark mass (the mass regularization). Our results are in agreement 
with those given in |28j . We have verified that our amplitude is both UV and IR finite and it's gauge 
invariant with respect to all the gluons. Pentagon and box amplitudes are separately gauge invariant 
with respect to the electro- weak vector bosons. As a consistency check we have also verified that the 
amplitude vanishes in the large quark-mass limit 



4 Numerical Results 

We have computed the amplitude numerically using the real polarization vector basis. There are 48 
polarized amplitudes for the case of jZg while 72 for the cases of ZZg and WWg. Given the number of 
diagrams, the number of polarization combination and the length of the amplitude, we have to run the 
code in parallel environment using a PVM implementation of the VEGAS algorithm |30|, 131] . We have 



4 



used more than 30 cores to run the code in parallel environment. Still it takes more than eight hours to 
get suitable cross-section for a given process. 

In such loop calculations, there is an issue of numerical instabilities due to large cancellations in the 
evaluation of the scalar and tensor integrals. To minimize this problem, we have used the OneLoop 
implementation of the scalar integrals, So, Co, and Dq. Even with this implementation, we have faced 
numerical instabilities in the evaluation of pentagon-tensor integrals Q This is related to inaccurate 
evaluation of Gramm determinants for certain phase-space points. These are the phase-space points 
where linear independency of external momenta (modulo 4-momentum conservation) does not hold to 
a very good accuracy. Indeed they lie on the kinematic boundaries of various channels. We judiciously 
throw away such points by making a gauge invariance test on the full amplitude. Since such points are 
few, given the nature of the integral evaluation, we don't expect significant impact of this aspect of the 
calculation. This is a common practice in such loop calculations. More complete discussion on this can 
be found in [321 S3]. 
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Figure 2: Variation of the cross-section with the Figure 3: Variation of the cross-section with the 
centre-of-mass energy for gg — > VV'g. centre-of-mass energy for gg — > VV' . 

In Fig. Q, we present the results of the cross-section calculation for three processes. These results 
include following kinematic cuts: 

P^' z ' w ' j > 30 GeV, \rf' z > w ' j \ < 2.5, R(j,j) > 0.6. 

We have also chosen factorization and renormalization scales as /i/ = /x/j = Pj,' Z ' W , as appropriate. 
We have used CTEQ6M parton distribution functions to obtain these results [53]. We note that at 
typical LHC energy, the cross-sections would be of the order of 100 fb. For example, at the centre- 
of-mass energy of the 8 TeV, the cross sections are 46.7 fb, 95.5 fb, and 225.2 fb respectively for the 
^Zg,ZZg and W + W~g production. Therefore, one may expect a few thousand of such events at the 
end of present run. But a W/Z boson is observed through its decay channels. If we consider the case 
when all the W/Z bosons are seen through their decays to the electron/muon only, then, with 20 fb -1 
integrated luminosity, the number of events for these processes will be approximately 60, 10, and 220 
respectively. However, if we allow one of the W/Z boson to decay hadronically, then the number of 

1 We would like to emphasize that we have also calculated di-vector boson production via gluon fusion for comparison 
as discussed below. These processes involve tensor integrals up to box type only and we have not faced any numerical 
instability in the calculation. The reduction method adopted in our calculation is indeed very stable numerically. 
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events would increase significantly. At the 14 TeV centre-of-mass energy, the numbers will be about a 
factor of three larger. The relative behavior of cross-sections of the three processes, as the center-of-mass 
energy varies, is quite similar to the case of di- vector boson production via gluon fusion as shown in 
Fig. g). This common behavior is mainly due to the relative couplings of the processes listed above in 
(3). We find that at 14 TeV the cross-sections for our processes are 20 — 30% of those for corresponding 
di-vector boson production (without jet) processes. We can also compare the contribution of the loop 
processes with the tree-level processes. We find that the processes gg — > jZg,W + W~ g make about 
4 — 5% contribution to the processes pp — > jZj,W + W~j, while gg — > ZZg makes a contribution of 
about 10 — 15% to the pp — > ZZj process. Here 'j' stands for a jet. This is quite similar to the case 
of di-vector boson production. Tree level estimates are obtained using MadGraph [35]. We have also 
compared our results for the W + W~j production with those of Melia et al. [20]. Though they have 
considered the leptonic decays of VT-bosons and kinematic cuts, choice of scales and parton distributions 
etc are quite different, the percentage contribution of gluon-gluon channel as compared to the LO cross- 
section is same within the allowed range of uncertainty, i. e. 4 — 5%. The values of the cross-section are 
more strongly dependent on the values of parameters and kinematic-cuts; still two results are similar if 
we take into account quoted uncertainties and branching ratios. The contribution of these gluon fusion 
processes can be even larger in appropriate kinematic regime. 

We now discuss few kinematic distributions at 8 TeV collider centre-of-mass energy. These distribu- 
tions remain same, characteristically, at 14 TeV centre-of-mass energy. The invariant mass distributions 
for pairs of vector boson are shown in Fig. Q. The positions of the peaks are related to the masses of 
the vector bosons in each case. In Fig. ^ the transverse momentum distributions for gluon is given for 
three processes. The major contribution to the cross-section comes from low Pt region, as we expect. 
The cross-section is very sensitive to Pt cut on gluon as it can come from one of the incoming gluons 
in the box contribution. We have compared the Pt distributions of 7 and g for the case of ^Zg in Fig. 
Q. The rapid fall of P T distribution compared to that of Pj. can be seen clearly. This is because the 
photon is coming from the quark loop, rather than from a quark- leg (as in the case of tree-level process). 
A larger cut on this transverse momentum will enhance the relative contribution of the loop process. 
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Figure 4: Invariant Mass distribution of vector boson Figure 5: Transverse momentum distribution of Glu- 
pairs at 8 TeV. ons at 8 TeV. 
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Figure 6: Comparison of distributions of 7 and g for "fZg production at 8 TeV. 



5 Conclusions 

We have presented the results of the cross-section calculations for the processes gg — > ^Zg, ZZg and W + W~g. 
These are one-loop processes and proceed via pentagon and box diagrams. We have kept only jZg like 
contributions for the cases of ZZg and WWg. There is decoupling of the top-quark. For the WW g 
case, the contribution from the first two generations of the quark is included. We have made a number 
of checks on our calculation. We have verified the cancellation of UV and mass singularities. We have 
also checked gauge invariance with respect to all the gauge particles. Axial-vector contribution does not 
interfere with the vector contribution and it is separately gauge invariant. We have observed qualitative 
similarity of these processes with the corresponding di- vector boson production cases. All these processes 
have already generated several hundred events at the LHC. In appropriate kinematic regime one should 
be able to observe these processes. The contribution of the process gg — > ZZg is specially significant. 
It may have bearing on the searches of the Higgs boson and new physics scenarios that lead to two Z- 
boson in the final state. Detailed and more complete study of these processes will be reported in future 
publications. 

AS would like to acknowledge fruitful discussions with R. Frederix and D. Zeppenfeld. He would also 
like to thank Dr. S. Veseli for his help with AMCI package. 
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